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ABSTRACT 

We use the galaxy cluster X-ray temperature distribution function to constrain the 
amplitude of the power spectrum of density inhomogeneities on the scale correspond- 
ing to clusters. We carry out the analysis for critical-density universes, for low-density 
universes with a cosmological constant included to restore spatial flatness and for gen- 
uinely open universes. That clusters with the same present temperature but different 
formation times have different virial masses is included. We model cluster mergers us- 
ing two completely different approaches, and show that the final results from each are 
extremely similar. We give careful consideration to the uncertainties involved, carrying 
out a Monte Carlo analysis to determine the cumulative errors. For critical density our 
result agrees with previous papers, but we believe that the result carries a larger un- 
certainty. For low-density universes, either flat or open, the required amplitude of the 
power spectrum increases as the density is decreased. If all the dark matter is taken to 
be cold, then the cluster abundance constraint remains compatible with both galaxy 
correlation data and the COBE measurement of microwave background anisotropics 
for any reasonable density. 

Key vifords: cosmology: theory - dark matter. 



1 INTRODUCTION 



1994 



Gorski et al. 1995; Liddle et al. 1996a), which possesses 



One of the most important constraints that a model of large- 
scale structure must pass is the ability to generate the cor- 
rect number density of galaxy clusters. This is a crucial con- 
straint, because a considerable amount is known about clus- 
ters from both optical and X-ray measurements. As clusters 
are relatively rare objects, in the standard picture they must 
form from fairly high peaks in the original density field and 
hence their abundance is very sensitive to the normalization 
of the power spectrum on those scales. As a result, the clus- 
ter abundance has been cited (White, Efstathiou & Frenk 
1993) as one of the strongest pieces of evidence against the 
standard cold dark matter (CDM) model when that model 
is normalized so as to reproduce the microwa ye background 



aniso t ropies as seen by the COBE satellite ( Bennett et al 
19961 ; [Banday et al. 1996| ; ^ orsk ct al. 199w Hinshaw et al 



1996) 



Amongst the many ways in which the CDM model can 
be salvaged, a popular choice has been to reduce the to- 
tal matter density. The most common context within which 
such models have been studied is in spatially flat cosmolo- 
gies, where a cosmological constant has been introduced to 



somewhat different late-time dynamics and also a different 
COBE normalization. 

In this paper, we re-assess the cluster abundance con- 
straint in CDM models, encompassing spatially flat models 
with both a critical and sub-critical matter density and also 
open universe models. We use the X-ray temperature distri- 
bution function, taking advantage of a number of extensive 
hydrodynamical cluster simulations which have now been 
carried out to calibrate this to Press-Schechter theory. 



2 THE POWER SPECTRUM 

Except for a brief discussion on COBE near the end of the 
paper, we shall be interested in scales short enough that 
even in the open case one can consider space to be flat. 
This allows us to specify the power spectrum of t he de nsity 
contrast, follo wing t he notation of Liddle & Lyth (1993) and 
Liddle et al. (1996a), as 



make up the 


shortfall 


in the matter density (Efstathiou, 


Sutherland & 


Maddox 


199C ; Kofman, Gnedin & Bahcall 


199J; Klypin, 


Primack 


& Holtzman 1995). However, re- 
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cently attention has also been focused on the possibility 
of a genuinely open cosmological model ( Ratra fc Peeble; 



where k is the wave number, a is the scale factor, H = a/a is 
the Hubble parameter, il is the density parameter, subscript 
'0' indicates the present value, T{k) is the transfer func- 
tion and the factor g{Q), defined later, accounts for the rate 
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of growth of density perturbations relative to the critical- 
density case whose growth is given by the (aH)'^ factor. The 
quantity 5^{k) specifies the shape of the primordial spec- 
trum; it is a constant for the Harrison-Zel'dovich spectrum, 
and 5h oc for a 'tilted' spectrum with spectral index n. 
A quantity related to the power spectrum is the dispersion 
a{R) of the density field smoothed on a (comoving) scale R, 
defined by 

dk 



a\R) 



W\kR)Vs{k)- 



k 



(2) 



To carry out the smoothing, we shall always use a top-hat 
window function W{kR) defined by 



(3) 



A large galaxy cluster has a mass of around 10^^ Mq; 
the amount of matter required to make such a cluster 
would originally (before collapse) be contained in a comov- 
ing sphere of radius 9.5(/i/no)^'^^/i~^ Mpc, where h is the 
present Hubble constant in units of lOOkms"^ Mpc~^. Tra- 
ditionally, the cluster abundance is used to place a constraint 
on the dispersion of the density contrast at the scale 8/i~^ 
Mpc, denoted erg, and some assumption regarding the shape 
of the power spectrum is required to shift from the scale 
actually constrained by the observations on to this one. For 
CDM models as considered here, the p ower spectrum is ac- 
curately given by Bardeen et al. (1986) as 

ln(H-2.34q) 



7cDM(g) 



2.34g 



(4) 



[l + 3.89g 4- (f 6.1g)' + (5.46g)^ + (6.71q)*] 



-1/4 



with q = k/hr, w here the 'shape parameter' F is defined as 
(Sugiyama 1995) 

(5) 



r — f2o/i exp(— SIb — Qb/^o) 



with f^B the baryon density. This fit to F is good for both flat 
and open universes. However, we have no particular need for 
it, as we can simply specify models by F. If the spectral in- 
dex n equals 1, then, provided that F is chosen in the range 
F = 0.230lo'g34 (these limits being 95 per cent confidence), 
a reasonable fit to the shape of the galaxy correlation func- 
tion is obtained. This allowed interval for F was calculated 



using the data points in table 1 of Peacock & Dodds ( 1994 ) , 
omitting the four corresponding to the smallest scales, as 
at these scales the working assumption of a linear bias used 
in their calculation seems to break down (Peacock 1996). 
Because galaxy correlations are presumed to be biased rel- 
ative to the mass, this fit tells us nothing about the overall 
normalization of the matter power spectrum. 

It is simplest to use an approximation to the true sh ape 
of a{R) in the vicinity of 8/i"^ Mpc. White et al. (|l993a|) use 



a power-law fit where the exponent is related to F. Because, 
particularly for low densities, we require accuracy over a 
greater range of scales, we use the more accurate fit 

/ R \ "^'^^ 

^(R) - ^T-i\j- , (6) 



8/1-1 Mpc 



where for the CDM spectra we adopt the form 



7(i?) = (0.3F-f-0.2) 



2.92 -f log 



R 



8ft- 1 Mpc 



(7) 



If one chooses a tilted primordial spectrum, n ^ 1, then 
the best-fitting F is changed, but the actual shape of the 
power spectrum is still obliged to fit the shape of the galaxy 
correlation function and so the results that we obtain will 
not depend significantly on n; indeed, even the F dependence 
will prove unimportant provided that it is restricted to lie 
within the favoured range. 

In a critical-density universe, clusters are anticipated 
to have formed very recently. As we shall see, in low-density 
universes clusters will have formed much earlier, and we shall 
need to take into account the redshift dependence of the 
power spectrum. In CDM universes, the shape of the power 
spectrum is unchanged at low redshift, so we need only con- 
sider the redshift dependence at a single scale which we take 
to be 8h~^ Mpc. This dependence is different depending on 
whether we have an open model (which we label 'OCDM') 
or a flat model ( which we label 'ACDM'). Following Carroll, 
Press & Turner (1992) we introduce a growth suppression 
factor (/(n), as used in equation (|l|). This gives the total 
suppression of growth of the dispersion relative to that of a 
critical-density universe, and is accurately parametrized by 



gm = 



1 + 
1 

70 



4/7 



+ 



140 140 



4/7 



(OCDM); (8) 



(ACDM). (9) 



These formulae can be applied at any value of fl. For a 
matter-dominated low-density universe, the redshift depen- 
dence of SI is given by 

1 + z 



1 + floz 



(1+^) 



1 - fio + (1 + z)^Qo 



(OCDM) ; (10) 
(ACDM) . (11) 



Since the growth law in a critical-density universe is erg (z) oc 
{1 + z)~^, the redshift dependence for arbitrary fio is there- 
fore 



erg (2) = (78(0) 



g{niz)) 1 



(12) 



gi^o) 1 + z' 

using the appropriate formulae for g{Q,) and ^l{z) depending 
on whether the universe is open or flat. 



3 PRESS-SCHECHTER THEORY 
3.1 Number densities 

The most accurate way of assessing the cluster abundance is 
via numerical simulation. However, there is an excelle nt an- 
alytic alternative which is the Press-Sc hechter theory (Press 
& Schechter 1974| 



Bond et al. 1991). For the case of a 



critical-density universe, this has been extensively tested 
against TV-body simu lations and found to fare extremely well 
(Lacey & Cole 1994). Less attention has been directed to- 



wards testing the theory in open universes, but those studies 
that do exist suggest that it continues to work well, as one 
might expect since its derivation relies only on statistical 
arguments. 

Press-Schechter theory states that the fraction of ma- 
terial residing in gravitationally bound systems larger than 
some given mass M is given by the fraction of space in which 
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the linearly evolved density field, smoothed on that mass 
scale, exceeds some threshold Sc. Originally the choice of Sc 
was motivated via the spherical collapse model, but it should 
now be calibrated via A''-body simulations. The appropriate 
formula is 



n{> M{R),z) 



• erfc 



V2aiR,z) 



(13) 



where 'erfc' is the complementary error function. Here R 
is the comoving linear scale associated with M, R^ — 
3M/47rpb, where pb is the comoving background density 
which is constant during matter domination. In equation 
(p^, the right hand side has been multiplied by a factor 
two to allow the material in underdense regions to partici- 
pate; in original treatments this seemed very ad hoc but it 



has now been more or less well justified (Peacock & Heavens 



199C ; Bond et al. 1991), and in any case it is incorporated in 
the TV-body calibration. The required value of Sc is quite de- 
pendent on the choice of smoothin g window used to obtain 
the dispersion (Lacey & Cole 1994). We shall use a top-hat 
window function. Were one to use a Gaussian window the 
threshold would be significantly lower (or alternatively the 
associated mass could be increased as discussed by Lacey & 
Cole 1994); claims in the literature of a large uncertainty 
in the Press-Schechter calculation are often due to quoting 
thresholds for different types of smoothing and one should 
always be careful to specify which is being used. 

For spherical collapse, the appropriate threshold is 5c = 
1.7±0.1. In the generic non-spherical situation, one must be 
careful to specify what is meant by collapse; if one considers 
collapse along all three axes the threshold is higher, whereas 
collapse along the first axis (pancake formation) or the first 
two axes (filament formation) corresponds to a lower thresh- 
old (Monaco 1995). Since large clusters are relatively rare, 
one can assume that shear did not play an important role 
and hence t heir collapse can b e considered to have occurred 
spherically (Bernardeau 1994). 

One might expect Sc to depend on the background cos- 
molog y. However, this seems not to be the case — Lilje 
( [1992|), L acey & Cole ( |l993| ) and Colafrancesco & Vitto- 
rio (1994) found that, at least for any type of collapse where 
(5c < 1.7 in a fiat universe, the value of Sc varies at most by 
5 per cent when one goes from a critical-density universe to 
one with Qq =0.1. Since at earlier epochs low-density uni- 
verses become closer to critical-density ones, the variation 
will be even less as one goes to higher redshifts. 

Equation ( ]l^ ) gives the total amount of material in 
structures above a given mass. However, we shall be in- 
terested in the number density of objects within a given 
mass range. The comoving number density of clusters in a 
mass interval dM about AI at a redshift z is obtained by 
differentiating equation ( p^ with respect to the mass and 
multiplying it by p\^/M. This gives 



n{M, z) AM = 

2 pb _ 

TT M a'' 



(14) 



Aa{R,z) 



{R,z) dM 



exp 



2a^{R, 



dM. 



Following White et al. ( 1993a ), we can simplify this for 
R in the vicinity of 8h~^ Mpc by using our analytic approx- 
imation to the power spectrum, equation IM), to calculate 



the derivative in equation ([l^). This gives 
n(M, z) dM = 



(15) 



2 Pb 2.92(0.3r + 0.2)(5c 
^JP 3a{R,z) 



exp 



5'c 



2a^{R, 



dM. 



3.2 Formation redshifts 



We shall be interested in the formation redshifts of the clus- 
ters that we see at the present epoch. The literature contains 
two very different ways of obtaining this (Lacey & Cole 1993, 
1994; Sasaki 1994), which each have their advantages and 
disadvantage s. In this paper we shall consider both. 

Sasaki (1994) showed through simple physical and 
mathematical arguments that the comoving number density 
of clusters with mass in a range dM about M, which virial- 
ize in an interval dz about some redshift z and survive until 
the present without merging with other systems, is given by 



N{M,z)dMdz ■ 



(16) 



(5c n{M,z) da{R,z) 
'cr^{R,z) a{R,z) di 



a{R,z) 
a{R,0) 



dM dz. 



The redshift independence of the shape of the CDM power 
spectrum allows us to write 



N{M, 



dMdz = 
S' 



(17) 



n{M,z) das{z) 



^-{R,z) aii{z) dz 



(78 (z) 

^8(0) 



dM dz , 



where (ys{z) and d(j%(z)/dz are calculated using equation 
(p^. In equation (^^, the expression within the square 
brackets gives the formation rate of clusters with mass M at 
redshift z, whereas the factor outside gives the probability 
of these clusters surviving until the present. The approxi- 
mation leading to this equation was to assume that the ef- 
ficiency of the destruction rate of clusters through mergers 
has no characteristic mass scale, so that merging proceeds 
in a self-similar way along the entire mass range. Though we 
do not expect this to happen for physically motivated power 
spectra, such as the ones under consideration, it should be 
a fairly good approximation across a limited range of scales. 
In order to be consistent with the Press-Schechter formal- 
ism, it then turns out that the efficiency of the destruction 
rate must be only a function of redshift. Blain ^ Longair 
(1993a,b) also worked within the Press-Schechter frame- 
work, instead assuming various physically reasonable merger 
probability distributions, and they found results numerically 
similar to Sasaki's. 

An alternative approach is that of Lacey & Cole (1993, 
1994), who attempted to construct a merging history for 
dark matter haloes based on the random walk trajectories 
technique. This approach is much closer to physical reality 
than that of Sasaki, but necessarily much more complicated. 
They considered two alternative possible techniques for com- 
puting the distribution of halo formation times, one based on 
an analytical counting method and the other based on Monte 
Carlo generated merging histories; a comparison they made 
with results from A'^-body simulations shows clearly that the 
first method is better, providing a good fit to the A-body 
results. In general, the calculation of the distribution of halo 
formation times using the analytical counting method has to 
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be performed numerically, but for a white noise power spec- 
trum, n = 0, an analytic expression is available. Usefully, 
it turns out that the distribution of halo formation times is 
almost independent of the shape of the power spectrum, so 
one can use the analytical expression for n = as a good 
approximation to the numerical calculation for any n. 

A drawback of the Lacey & Cole approach is that one 
has to choose when to consider a given dark matter halo to 
have formed, since the last infinitesimal amount of mass is 
still being accreted at the present. That is, at what fraction 
of its final mass is a halo to be considered to have effectively 
formed? Bearing in mind that difi'erent properties of a clus- 
ter may be affected to different degrees by increasing the 
cluster mass, clearly the definition of when a cluster has ef- 
fectively formed depends on the particular cluster property 
under study. A cluster will have effectively formed if the 
property barely changes until the cluster reaches its final 
mass. According to Lacey & Cole (1993, 1994), the proba- 
bility that a cluster with present virial mass M would have 
formed at some redshift z is then given by 



matter density. Bartlett & Silk (1993) tested a variety of 



where 



dw{z) 
dz 



(18) 



p(w{z)) = 2w{z) - 1) erfc (^^^ - (19) 



(/ - 2) exp ' 

TV ^ ' 



and 



10(2) = 



4(a(M,0)/g(M,z)^ 1) 
v/a2(/M,0)-a2(M,0) 



(20) 



with / the fraction of the cluatcr maaa aaacmblcd by rcdahift 
As t he shape of the CDM power spectrum is redshift in - 



dependent we have a(M, 0)/(j(M, z) = as{0)/o^z), where 
as{z) is given by equation (y_2|). Expression ( |l9| ) was ob- 
tained for a power-law spectrum with index n — 0, but, as 
mentioned before, numerical results show that p{w{z)) de- 
pends only very weakly on n (Lacey & Cole 1993) so we will 
use it even though, at the scales we are interested in, the 
spectral index is closer to n = —2. 



4 THE CLUSTER ABUNDANCE 

The first attempts to constrain the power s pectr um via the 
cluster abundance w ere m ade by Evrard (1989) and then 



by Henry & Arnaud (1991). They both considered only the 



critical-density case, Evrard using the velocity dispersion to 
determine the mass via the virial theorem while Henry & Ar- 
naud used the X-ray temperature. Both obtained very simi- 
lar results, though the quoted errors of the latter were much 
smaller . Oth er pie-COBE analyses were made by B ond fc 
Myers ( |l992| ), by Bahcall & Cen (1992, 1993), by Lilje ( |l992| ) 
who consid ered low-density fiat models and by Oukbir & 
Blanchard (1992) who discussed the open case. Subsequent 



to the COBE observations. White et al. ( 1993a ) carried out 



an analysis with extra input from A^-bod y simu lations, ob- 
taining a si milar result again to Evrard (1989) and Henry 
& Arnaud (1991) in the case of critical density and extend- 
ing the calculation to the case of flat universes with a low 



flat-space models against the data, using the then current 
flrst year COBE normalizatio n which is some w ay below 
that currently recommended ( |G6rski et al. 1996 ). Balland 
& Blanchard (1995) considered the hot dark matter case. 
Recently, Liddle et al. (1996a) briefly described a new cal- 
culation in the case of an open universe. That calculation 
extended the type of analysis made earlier, by attempting 
to take into account that clusters with equal mass which 
virialize at different redshifts have distinct properties, such 
as velocity dispersion and X-ray temperature, at the present. 
As well as providing a more detailed account of that open 
universe calculation, in this paper we carry out a similar 
analysis for flat universes. 

When applied to rare objects, the number density pre- 
dicted by the Press-Schechter theory is extremely sensitive 
to the dispersion a{R). This is of great advantage, because 
it means that even if the number density is not well known 
the error this induces in estimating as is small. Much more 
crucial in this application is to establish estimates of the 
cluster masses as accurately as possible. There are presently 
three methods in use for mass estimation. The cluster ve- 
locity dispersion provides one such means. Unfortunately, 
cluster catalogues assembled from optical data are prone 
to contamination from foreground and background galaxies 
mis-identified as part of the cluster, and may also be affected 
by projection effects and by the possibility of velocity bias. 
An alternative means of cluster identification is via X-ray 
emission from the gas resting in a deep potential well. Since 
X-rays are produced only in clusters, there are no problems 
with foreground and background contamination (unless two 
clusters happen to lie on top of one another in projection, 
and even then discrimination may be possible as luminos- 
ity goes as a steep power of the cluster richness (Henry & 



Arnaud 1991)). The final method is via weak lensing of back- 
ground galaxies ( Kaiser fc Squires 1993| ) . This promises to be 



a very interesting technique for the future, but at present has 
not been widely applied. Consequently, we choose to adopt 
the X-ray data. 

The observed number density of clusters per unit tem- 
perature, n(fcBT), at z — has bee n det ermined by Edge et 
al. ( [l990| ) and by Henry & Arnaud ( |l99l| ). These are in good 
agreement; we shall use the latter^ We shall concentrate on 
large clusters by considering those with mean X-ray temper- 
ature 7 keV; those are observed to have a present number 
density per unit temperature interval of 



n(7 keV,0) = 2.0^1°^ x 10''' h'^ Mpc'^ keV"^ , 



(21) 



In order to apply the Press-Schechter formalism, one 
needs to relate the X-ray temperature to the virial mass. As 
we shall see, in the case of a critical-density universe one 
expects that all clusters formed fairly recently and there is 
more or less a one-to-one correspondence between a given 
temperature and a given virial mass. In low-density uni- 
verses, clusters can form much earlier, and the clusters we 
see today of a given temperature would have formed at a 
range of different redshifts. Lilje ( 1992| ) (see also Hanami 



* Eke, Cole & Frenk (1996) have recently pointed out two errors 
in the derivation by Henry & Arnaud, which fortunately almost 
exactly cancel each other out. 
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1993) has demonstrated that clusters of the same present 
temperature, but different formation times, will in general 
have different virial masses in accordance with the scaling 
relation 



where 



3/2 



(22) 



-3/2 



+ 



-3/2 



(fceT) 



3/2, 



X = 



(23) 
(24) 



with ^ the ratio between the cluster and background densi- 
ties at turnaround, and 



2 + e - 77/2 ' 



(25) 



where Tm and are the radii of turnaround and virialization 
respectively. The parameter e represents the difference in the 
total energy of a cluster from the one obtained by assuming 
the cluster to be an ideal virialized system collapsed from 
a top-hat perturbation. The two main processes by which 
such a difference can be introduced are radiative cooling and 
dynamical relaxation due to the presence of substructure. In 
both cases e > and there is a loss of energy by the cluster, 
thus making it more compact. Whilst the first is generally 
regarded as unimportant due to the fact that the cooling 
time of a typical rich cluster is larger then the age of the 
universe, the second could have a significant impact on the 
cluster final state. However, because the discussion on the 
importance of this process is still wide open, we choose to 
use e = 0. In the above expressions Ao = \/3Hq, so Ao = 
1 — fio in a flat universe. Also, Zc and Zm are respectively 
the redshifts of cluster collapse and turnaround; they are 
related by the fact that the time of collapse and virialization 
tc is twice the time of turnaround tm, as the expansion and 
subsequent collapse of a cluster are symmetric about the 
time of turnaround. We shall need to calculate Z m given 
Zc- The redshift -time relation for open models is (Kolb & 
Turner |l990|) 



^0 



° 2(1-^0)3/2 

2(l-no)^/^(noz+l)i/^ 
no{l + z) 



(26) 



cosh 



- fio + 2 

z + Qo 



.1 /S2o£_ 

V ^0 



and for flat models is (Charlton & Turner 1987) 



In 



3A, 



/2 • 



(27) 



[,\o(i + .)-=']^/" + [Ao(i + .)--^ + r i o]^/" 



1/2 



In the limiting case of a critical-density model we have 

t = lH^\i + zr'^\ 

As we know that tm = tc/2 we can therefore calculate Zn 
from the above equations through an iterative procedure. 



(28) 



The parameter ^ can be obtained by solving the equa- 
tion of motion of the outermost mass shell of a cluster, and 
using the fact that at turna round this shell has zero veloc- 
ity. Following Hanami ( 1993| ) (see also Martel 1991) we then 
have 



where 



(1 + Zm) 



-3/2 



1/2 



Ao 
Ao 



i-(i + C)y + Cy« 



Ay (OCDM); (29) 

1/2 

Ay (ACDM),(30) 



(31) 



Whilst for open models the integral in equation (|29|) can be 
solved analytically, yielding 



r(l + . 



(32) 



for fiat models the integral in equation ( ^ ) has to be solved 
numerically. Hanami (199J) provided a single fitting func- 
tion which proves good for the open case but not nearly as 
good for the flat case. We calculate improved fitting func- 
tions by using the trick that any epoch can be regarded 
as the present provided one uses the appropriate value of 
Q,{z) as given by equation (^o|) or (|l^). Taking the epoch 
of turnaround to be the present leaves x depending only on 
the value of Q,, and we can then fit a simple function to the 
numerical result. We find 



-tm 



where Q = Q{zm), and within an error of 2 per cent 



fm = 0.76- 
/(«) = 0.73- 



0.25 
0.23 



(OCDM) ; 
(ACDM) . 



(33) 

(34) 
(35) 



The scalings in equation ( p2| ) have been tested by hy- 
drodynamical N-hoAy simulations in the case of a critical- 
density model where they have been found to hold very well 
(Navarro, Frenk & White 1995). 

The crucial question is the value of the proportionality 
constant in equation (|2^). We fix this by taking advantage 
of a set of hydrodynamical A'^-body simulatio ns carr ied out 
in the critical-density case by White et al. ( [l993bD . They 
created a catalogue of twelve simulated clusters with dif- 
ferent temperatures, and found that a cluster with present 
mean X-ray temperature 7.5 keV corresponds to a mass 
Ma = (1.10 ± 0.22) X 10^^ h'^ Mq within the AbeU ra- 
dius (1.5/i~^ Mpc) of the cluster centre. The quoted error is 
the la dispersion within the catalogue. The conversion from 
an AbeU radius to the virial radius is then standard, via the 
result that the simulated clusters have a density profile in 
the outer regions described by Pc(r) tx )-~2'4±o i ( yy^te et 



al. 1993b 



Metzler & Evrard 1994; Navarro et al. 1995). In 



a critical-density universe the virial radius encloses a den- 
sity 178 times the background density, and hence through a 
Monte Carlo procedure, where we assume the errors in Ma 
and in the exponent of Pc{r) to be normally distributed, we 
obtain Mv = (1.23±0.32) x 10^^ Mq for the virial mass 
corresponding to a 7.5 keV cluster. The remaining uncer- 
tainty is th e virialization redshift o f such a cluster, whic h is 
estimated ( Metzler fc Evrard 1994 ; Navarro et al. 1995 ) as 
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Zc — 0.05 ±0.05. This enables the normahzation of equation 
(p^, and for the particular case of a 7 keV cluster, the ob- 
servation that we are using, one then has a virial mass given 
by 



(1.2 ±0.3) X 10'" ^0 

3 



1/2 



3/2 



n 



X 

3/2 



-1/2 f2^y X (36) 

(l±2c)"^^'fe"' Mq . 



We shall now compute the present comoving number 
density of clusters for a unit temperature interval about 7 
keV in a given cosmology, in order to compare it with the 
observed value. We will do t his in two different ways, one 
using the results of Sasaki (1994) and the other those of 
Lacey & Cole (1993, 1994). 



4.1 The Sasaki method 



According to Sasaki (1994), the comoving number density 



of clusters that virialize at redshift z with some virial mass 
M, and that then survive up to the present without merging 
with other systems, is given by equation (|l^) . For our appli- 
cation, we wish to consider at each redshift the appropriate 
virial mass which gives rise to a 7 keV cluster, via equation 
( p6| ) with Zc = z. Using the chain rule, we can obtain the 
comoving number density of clusters per temperature inter- 
val d(fcBT) that virialize at each redshift z and survive up 
to the present such that they have a present mean X-ray 
temperature of 7 keV: 

'^^'^ -N{M,z)d{kBT)dz 



N{kBT,z)d{kBT)dz = 



d(fcBT) 
3 M 



2 ksT 



N{M,z)d(kBT)dz, (37) 



where the second equality uses equation 
the final expression 



22h. This yields 



N{kBT,z)d{kBT)dz = 
3 M 



(38) 



n{M,z) das{z) 



2 fcBTcr2(i?,z) (78(2 = 0) dz 



d(fcBT) dz . 



In order to compute the present comoving number density 
of clusters per unit temperature at 7 keV, one needs to in- 
tegrate this expression from redshift zero up to infinity with 
the appropriate cosmological information inserted. Fig. 1 
shows plots of the integrand as a function of redshift, thus in- 
dicating when the surviving clusters predominantly formed, 
for a critical-density model and for both a flat and an open 
model with Qo ~ 0.3. From the integrated expression, one 
can find the required erg to obtain the observed number den- 
sity. 

For both open and flat models, one can quote the result 
in the form of the required as as a function of flo- The best- 
fitting value is given by 



0-8 



0.60 ^o'^*"''^ 



where 



0.34 ±0.28 Ho -0.13 
0.57- 0.13 ±0.08 0,1 



(39) 

(OCDM) , (40) 
(ACDM) (41) 



Q„=1.0, CDM 
n;,=0.3, OCDM 
Q„=0.3, ACDM 




are fitting functions representing the changing power-law in- 



Figure 1. This shows the differential rate of formation of clusters 
with mean X-ray temperature 7 keV that survive to the present 
using the Sasaki method, where we take d(fcBT) to be the unit in- 
terval. Curves are shown for three cosmological models: a critical- 
density universe (solid), an open universe with Qo = 0.3 (dashed) 
and a flat model with Qo = 0.3 (dotted). Each curve is normalized 
to produce the same present number density. 



dex of the Qo dependence and are accurate within 2 per 
centfj. 

While the central value is itself of interest, it is vital 
to know what range of values of as about this is permitted. 
The uncertainties arise from a variety of sources. Ranking 
them in order of importance as contributions to the overall 
uncertainty in as, we have first the uncertainties in relat- 
ing the cluster temperature to its virial mass in equation 
(M) and the uncertainty in the Press-Schechter threshold 
parameter 5c, then the uncertainty in the observed num- 
ber density in equation (^), and finally the uncertainty in 
the observational value of T. Since the analysis is both nu- 
merical and non-linear, we estimate the errors via a Monte 
Carlo procedure, whereby we model the different uncertain- 
ties as Gaussians (the observed number density is modelled 
as a Gaussian in its logarithm, as suggested by the error 
bars) . Realizations are then drawn from the Gaussian distri- 
butions and processed through the calculation to determine 
the required as; the distribution of these is then taken as 
the uncertainty in erg. While this procedure implicitly imag- 
ines the errors to be statistical whereas in reality they may 
be predominantly systematic, the fact that there are sev- 
eral sources of errors, none of which dominates completely, 
means that this procedure should not be excessively strin- 
gent, and indeed our estimate of the uncertainty turns out 
to be rather larger than others that have appeared in the 
literature thus far. 

The uncertainty in the observed slope of the power spec- 
trum r is included in the overall uncertainty. We find that, 
for each fio between 0.1 and 1.0, the distribution of as is 
close to lognormal. In the OCDM case, the 95 per cent con- 
fidence limits are ±32 per cent and —24 per cent, indepen- 
dent of Slo to a good approximation. In the ACDM case, the 
uncertainty becomes larger at low density; a satisfactory fit 



t Note that, although we have made several alt eration s as com- 
pared to the OCDM calculation in Liddle et al. ( |l996£j ), the final 
result is very similar. 
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to this increase in uncertainty is to multiply the OCDM 
uncertainty by a factor f2°-^*'°sio'^o p^j. example, in a flat 
universe with Slo = 0.3 the 95 per cent confidence limits 
are increased to +37 per cent and —28 per cent, and the 
uncertainty climbs rapidly if fio is further reduced^ 



4.2 The Lacey & Cole method 

We will now use the results from Lacey & Cole (1993) to 
obtain an alternative estimate of as{Q,o) for open and fiat 
cosmological models. According to equation (llq), one can 
estimate the fraction of present clusters with virial mass M 
that formed at a redshift z. As we want to count only those 
present clusters with a mean X-ray temperature of 7 keV, 
for each M the corresponding value of z is uniquely fixed by 
equation (|36|). Therefore the product of the present comov- 
ing number density of clusters per unit mass with virial mass 
M (given by equation ( p^ ) with z = 0) with the fraction of 
those clusters with a present mean X-ray temperature of 7 
keV (given by equation ( p^ with z obtained from equation 
( p6| ) for the given M) uniquely defines the comoving number 
density of clusters per unit mass with present mean X-ray 
temperature 7 keV that formed at redshift z: 



N(M, z)dMdz = n(M, 2 = 0) p{z) dMdz . 



(42) 



Again using the chain rule, the present comoving number 
density of clusters per temperature interval d{kBT) with a 
mean X-ray temperature of 7 keV that formed at each red- 
shift z is then given by 



N(kBT,z)d{kBT)dz 



3 M 
2k^ 



n{M, 2 = 0) p{z) dMdz . (43) 



Once more the present comoving number density of clusters 
per unit temperature at 7 keV is obtained by integrating 
this expression from redshift zero up to infinity. 

As we mentioned in the last section, the approach of 
Lacey & Cole (1993, 1994) has the drawback of requiring 
one to define a criterion for when a cluster is effectively 
formed in terms of the fraction of the cluster final mass as- 
sembled at that moment, /. In their work Lacey & Cole 
defined the formation time of a cluster as the moment when 
half of its final mass has been assembled. For our purposes 
we require the moment after which any cluster mass increase 
leads to only a small change in the cluster temperature. To 
our knowledge the only comparison between the evolution of 
a clust er's m ass and its X-ray temperature is that of Navarro 
et al. (1995). We have already used their results to estimate 



the formation redshift of clusters in a critical-density uni- 
verse in order to normalize equation (|2^), where we took 
Zc ~ 0.05 ± 0.05. We can read from Navarro et al. (199. 



how much mass the clusters had assembled within that ap- 
proximate redshift interval. This gives / = 0.75±0.15, where 



t We also made an analysis where Sc was shifted to its 95 per 
cent limit and then the remaining errors analysed via the Monte 
Carlo procedure; this corresponds to treating the Sc uncertainty 
as entirely systematic. If one does this, the error bars at 95 per 
cent confidence are increased to -|-42 per cent and —31 per cent 
in the OCDM case, the relative increase being the same in the 
ACDM case. 




Q„=1.0, CDM 
n;,=0.3, OCDM 
Q„=0.3, ACDM 



Figure 2. This shows the differential rate of formation of clus- 
ters with present mean X-ray temperature 7 keV using the Lacey 
& Cole method, where we take d(fcBT) to be the unit interval. 
Curves are shown for three cosmological models: a critical-density 
universe (solid), an open universe with f2o = 0.3 (dashed) and a 
fiat model with f2o = 0.3 (dotted). Each curve is normalized to 
produce the same present number density. 



the error is supposed to correspond to a 2a confidence inter- 
val. This seems a physically reasonable result, and we shall 
assume that it remains true in all cosmologies. 

In Fig. 2 we show plots of expression (US) as a function 
of redshift, thus indicating when the clusters of temperature 
7 keV predominantly formed, defined as when 75 per cent 
of their mass had assembled. 

Integrating expression ( ^ ) from redshift zero to infinity 
and comparing the result with the observed number density 
( pl[ ) we once more obtain ag as a function of flo for both 
open and fiat models. Using the Lacey & Cole approach the 
best-fitting value is given with an error of less than 2 per 
cent by 

-c{no) 



0-8 = 0.60 fio 
with 

Ci^o) = 0.39-^0.23 00-0.23^0 
Ci^o) = 0.61 - 0.19 J7o 0.02^0 



(44) 

(OCDM) ; (45) 
(ACDM) . (46) 



The overall error in t he val ue of erg obtained using the 
method of Lacey & Cole (199J) comes from the sa me un - 
certainties that we had using the method by Sasaki (1994), 
plus the uncertainty in the value of the fraction of the cluster 
mass assembled at the time of effective formation. We again 
estimate the overall error via a Monte Carlo procedure, find- 
ing that the distribution of ag is still close to lognormal. In 
both the OCDM and ACDM cas es the results obtained us- 
ing the method of Lacey & Cole ( 199^ ) show an increase in 
the ove rall er ror in the value of erg relative to the method of 
Sasaki (1994). This increase becomes larger at low densities 
and is more important in the OCDM case. The 95 per cent 
confidence limits for the valu e of ag calculated using the 
method of Lacey & Cole (1993) are given to a good approxi- 



mation by +32D.l-"^°'^^'>^'> per cent and ~2'in[ 
cent in the OCDM case, and -1-3201] 



0.261ogior!o 



o.iTiogiono 





per 

per cent and 



-24no^"°*^io"'' per cent in the ACDM case. For example, in 
an open universe with Qo = 0.3 the 95 per cent confidence 
limits are -1-36 per cent and —27 per cent, and in a flat uni- 
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Figure 3. The dependence of eg on is shown as obtained 
using the two different methods of Sasaki (solid lines) and Lacey 
&; Cole (dashed lines). The lower two curves are for the open case, 
and the upper two for the flat case. 



verse, also with Ho = 0.3, the 95 per cent confidence limits 
are +38 per cent and —28 per cent. 



5 DISCUSSION 

In Fig. 3 we compare the values ob taine d for erg (Ho) using 
the two different methods of Sasaki (1994) and Lacey & Cole 



(1993). For OCDM the difference is less than 5 per cent down 
to n = 0.2, increasing to 11 per cent for = 0.1, whilst 
for ACDM the difference is less than 3 per cent down to 
Q. — 0.2, increasing to 9 per cent for O = 0.1. The difference 
between the results is much less than the other uncertainties. 
In order to present a definite result, we fit to the mean of 
the two methods to obtain our final result: 



ag = 0.60 
with 

C{no) = 



-C(no) 



0.36 + 0.31 no 
0.59- 0.16 



0.28 fit 
0.06 0,1 



(47) 

(OCDM) ; (48) 
(ACDM) . (49) 



Since the difference between the two methods is small, 
the overall 95 per cent confidence limits remain as before. 
The most conservative assumption is to c onside r those ob- 
tained using the method of Lacey & Cole (1993), which are 



given by +32f7°-'"°'^i°"° 



per cent and —240, 



TTiTlogiono 



per 



cent in the OCDM case, and +32f2Q 



0.261ogi(,no 



-24n 



,0.261ogiono 




per cent and 



per cent in the ACDM case. 



We do not intend a detailed comparison of these re- 
sults with other types of observation here (see Liddle et al. 
1995, 1996a, 1996b), but it is worth comparing the results 
with t he normalization from the four-year CO BE observa- 



tions 



([Bennett et al. 1996 ; Banday et al. 1996; [Corski et 



al. 1996; Hinshaw et al. 199(: ) for the case of scale- invariant 



(n = 1) initial density perturbations. These fix the ampli- 
tude of perturbations on large scales, 5^ as a function of 
0,0, independently of h to an excellent approximation. Fit- 
ting functions have been calculated by Bunn & Wh ite (i n 
preparation) using techniques from White & Bunn (1995), 
and are 




Figure 4. The cluster constraint (solid) is compared with the 
COBE normalization (dashed) for the range of F satisfying the 
galaxy correlation data. The upper figure shows the open case, 
and the lower one the flat case. In each, it is assumed that all 
the dark matter is cold. The upper COBE line corresponds to 
r = 0.27 and the lower one to T = 0.20. Both cluster and COBE 
constraints are at 95 per cent confidence. The central values are 
not shown. 



lo^^^K^o) 



3.76 a 



-1.57-0.10 Inflo 



(ACDM) . 



(51) 



IQ^^'SliOo) 



3.80 n, 



-0.70-0.38 In no 



(OCDM) : 



(50) 



The fitting functions are accurate within about 3 per cent 
for 0.2 < Oq < 1.0. The COBE normalization error bar on 
5h is 15 per cent at 2o. 

Although these values are independent of h, and indeed 
of the nature of the dark matter, they pick up a dependence 
on these when one computes the equivalent erg . We shall as- 
sume that the dark matter is all cold, so that we can use 
the r parametrization of the transfer function. The uncer- 
tainty in the observed value of F given after equation (^ 
propagates to the calculated value for ag when one uses the 
appropriate COBE normalization for some value of Q.o- We 
will therefore add in quadrature the uncertainty thus arising 
in erg and the 15 per cent error at 2a which appears in erg 
due to the COBE normalization uncertainty. 

The comparison is plotted in Fig. 4. The allowed range 
of values for Q.o goes down to about 0.20 in the OCDM case. 
In the ACDM case the cluster constraint remains compat- 
ible with COBE and the galaxy correlation function shape 
down to very low fio. However, in this latter case there are 
other reasons for believing that low values are not favoured, 
for example, values of fio below 0.4 require an anti-bias for 
optically selected galaxies due to their very high COBE nor- 
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malization. This is difficult to understand physically. The 
anti-bias of IRAS galaxies would be even more pronounced. 
However, this conclusion is driven by COB-B rather than the 
cluster abundance and can perhaps be partly alleviated by 
tilting the primordial spectrum. 

For completeness, in Fig. 5 we show the comoving num- 
ber density per unit temperature interval of clusters with 
mean X-ray temperature 7 keV that we should expect to ob- 
serve at redshift z, for three cosmological models normalized 
to the present central value. The Sasaki and Lacey & Cole 
methods give very similar results. Note that these are clus- 
ters which have that temperature at the given redshift; their 
X-rays would be redshifted on the way to us, hence making 
the apparent temperature smaller by a factor of (1 -I- z). 
As one progresses to a lower redshift, some of the clusters 
will merge and some new ones will form. At present the 
only complete surveys of clusters at high redshifts are flux- 
limited (e.g. Henry, Jiao & Gioia 1994). This poses various 
problems, the most serious of which is that as the cluster 
selection is then a function of their luminosity, not of their 
temperature, we need to know the temperature-luminosity 
relationship at the appropriate redshift to be able to com- 
pare the observations with our results. However, this rela- 
tion is very sensitive to several parameters like the cluster 
baryon fraction and the extent to which energy is injected 
into or removed from the intracluster medium, examples of 
the former being gas stripping from galaxies infalling into 
the cluster or supernova explosions in cluster galaxies, and 
of the latter being cooling flows. 

To recap on our main results, we have used a new 
method to compute the power spectrum constraint arising 
from the observed number density of clusters with a given 
X-ray temperature^. We have done this for a critical-density 
universe and for both flat and open low-density universes. 
Although we have assumed that all the dark matter is cold , 
the constraint would not change much ( Liddle et al. 1996h ) 
if one went, for example, to a cold plus hot dark matter 
model. In cases where such calculations have been done pre- 
viously, we support the previous results but typically find 
larger accumulated uncertainties. 

For the case of pure cold dark matter, we have also com- 
pared our results with constraints from COBE and from the 
galaxy correlation function (which constrains the shape pa- 
rameter r). We found that the cluster constraint is compat- 
ible with these in both the flat and open cases for any rea- 
sonable value of the density, failing only in the open case for 
Qo < 0.20. However, near critical density the corresponding 
h will be very low (below 0.5), and at the lowest permitted 
densities it will be very high (above 1.0). Further, in the 
spatially flat case the concordance of these constraints at 
low values of Qq will require optically selected galaxies to be 
anti-biased with respect to the dark matter distribution. 




•1.0, CDM 
•0.3, OCDM 
•0.3, ACDM 




Figure 5. The expected comoving number density of clusters for 
a unit temperature interval about 7 keV at redshift z for three 
cosmological models: a critical-density universe (solid), an open 
universe with CIq = 0.3 (dashed) and a flat model with Oq = 0.3 
(dotted). Each curve is normalized to produce the same present 
number density. The upper panel shows the result using the Sasaki 
merger calculation, and the lower one that of Lacey & Cole. 
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